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We present the first observations of large amplitude waves 
in a well-defined electron diffusion region at the sub-solar 
magnetopause using data from one THEMIS satellite. These 
waves identified as whistler mode waves, electrostatic soli- 
tary waves, lower hybrid waves and electrostatic electron 
cyclotron waves, are observed in the same 12-sec waveform 
capture and in association with signatures of active mag- 
netic reconnection. The large amplitude waves in the elec- 
tron diffusion region are coincident with abrupt increases in 
electron parallel temperature suggesting strong wave heat- 
ing. The whistler mode waves which are at the electron scale 
and enable us to probe electron dynamics in the diffusion re- 
gion were analyzed in detail. The energetic electrons (~30 
keV) within the electron diffusion region have anisotropic 
distributions with T e ±/T e \\ > 1 that may provide the free en- 
ergy for the whistler mode waves. The energetic anisotropic 
electrons may be produced during the reconnection process. 
The whistler mode waves propagate away from the center of 
the 'X-line' along magnetic field lines, suggesting that the 
electron diffusion region is a possible source region of the 
whistler mode waves. 



1. Introduction 

Magnetic reconnection is considered to be an important 
energy conversion process [Dungey, 1961] that occurs in a 
variety of plasma environments. At the Earth's magne- 
topause, it facilitates the entry of solar wind plasma and 
electromagnetic energy into the magnetosphere. Reconnec- 
tion sites are regions of strong wave activity covering a broad 
range of frequencies. Wave modes frequently observed near 
reconnection sites include the whistler mode (WH) waves 
[Deng and Matsumoto, 2001; Petkaki et ai, 2006], elec- 
trostatic solitary waves (ESWs) [Farrell et al., 2002; Mat- 
sumoto et al, 2003], lower hybrid (LH) waves [Cattell and 
Mozer, 1986; Bale et ai, 2002], kinetic Alfven waves [Chas- 
ton et al., 2005] and Langmuir/upper hybrid waves [Farrell 
et al, 2002]. The effect of different wave modes on the recon- 
nection process has been a problem of longstanding interest 
- for their role in anomalous resistivity, particle acceleration, 
energy transport and formation of reconnection sites [Huba 
et al., 1977; Labelle and Treumann, 1988; Treumann et al., 
1991; Drake et al, 2003]. 
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Observations of the electron diffusion region (EDR) have 
been made by Polar at the subsolar magnetopause [Mozer 
et al., 2002], by Wind in the magnetotail [Oieroset et al., 
2002] and by Cluster in the magnetosheath [Phan et al., 
2007] . Recent simulations and observations of EDRs during 
collisionless antiparallel reconnection in Earth's magnetotail 
[Ng et al, 2011] report that the diffusion region is charac- 
terized by a narrow extended layer containing electron jets. 
It is shown that the jets in the layer are driven by elec- 
tron pressure anisotropy P e || >> P e ± and the anisotropy is 
responsible for the structure of the EDR [Ng et ai, 2011]. 
Mozer [2005] has identified EDRs on the basis of the non- 
zero parallel electric field, a large perpendicular electric field 
compared to the reconnection electric field, a large electro- 
magnetic energy conversion rate and accelerated electrons, 
and a topological boundary that separates regions having 
different E x B/_B 2 flows with thickness of the order of the 
electron skin depth. In this paper, we concentrate on a spe- 
cific way to identify the EDR described by Scudder et al. 
[2012] who report spatially resolved diagnostic signatures of 
a demagnetized EDR observed by Polar at the Earth's mag- 
netopause. The five dimensionless scalar diagnostics that 
were used to find the EDR are peak electron thermal Mach 
numbers M e ± = < |y^ > > 1.5 where U e represents electron 
bulk velocity and < w e ± > is the electron thermal speed de- 
rived from the average perpendicular temperatures, electron 
temperature anisotropy An e = <T "^ > > 7, calibrated agy- 
rotropy of electron pressure tensor A(p e = 2 |*~"j > 1 where 
a = P e ±,i/P e ±,2, expansion parameters of guiding center 
theory indicative of demagnetization and strong (150eV) 
increases in electron temperature [Scudder et al, 2012]. 

This paper focuses on the WH waves which are an im- 
portant candidate for the anomalous resistivity, particle ac- 
celeration and heating. WH waves may be driven unstable 
by superthermal electrons with temperature anisotropics of 
T e x/T e || > 1 in the magnetosphere [Kennel and Petschek, 
1966] and current-driven plasma instabilities [Gurnett et al., 
1976] or energetic electron beams [Zhang et al, 1999] in the 
magnetotail. WH waves are one of the most ubiquitous 
wave modes observed in space plasmas. Observations of 
WH waves at the Earth's magnetopause have been made by 
Deng and Matsumoto [2001]. Electron anisotropy, due to 
compression of the magnetopause or LH drift waves, may 
be the generation mechanism of WH waves in the magne- 
topause current sheet [Karimabadi et al, 2004]. WH waves 
in the EDR may play a significant role in the microphysics 
of reconnection as they are excited on electron scales. It is 
believed that WH waves in the magnetopause current sheet 
may affect the instability of the current sheet to reconnec- 
tion via tearing. The generation of the out-of-plane compo- 
nent of the magnetic field is suggested to be a signature of 
whistler mediated reconnection [Mandt et al., 1994]. It has 
also been suggested that the strongest whistler emissions are 
observed on the most recently opened magnetospheric flux 
tubes due to magnetic reconnection [Vaivads et al, 2007]. 
One recent simulation study concludes that WH waves do 
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not control the dissipation processes of reconnection but are 
generated as a result of the reconnection processes [Fujimoto 
and Sydora, 2008]. In this paper, we present an example of a 
reconnection event at the sub-solar magnetopause observed 
by THEMIS. In section 2, we describe the data sets and 
analysis techniques. In section 3, we show the observations. 
Finally, we discuss the conclusions of our study in section 4. 



2. Data Sets and Analysis 

The THEMIS mission consists of five identically- 
instrumented spacecraft [Angelopoulos, 2008]. The Elec- 
tric Field Instrument (EFI) measures three components of 
the electric field [Bonnell et al, 2008]. The instrument 
provides continuous coverage at 128 samples/s in survey 
mode and waveform captures at 8192 samples/s in wave 
burst mode. The Magnetic Fields Experiment (MFE) in- 
cludes a Flux Gate Magnetometer (FGM) which measures 
DC magnetic field with a sampling rate of 128 samples/s 
in the high rate mode or 4 samples/s in the low rate mode 
[Auster et al, 2008] and a Search Coil Magnetometer (SCM) 
which measures magnetic fluctuations sampled at 8192 sam- 
ples/s in the burst mode [Roux et al, 2008; Le Contel et al., 
2008] . Particle data are measured by the Electrostatic Ana- 
lyzer (ESA) [McFadden et al, 2008] and by the Solid State 
Telescope (SST) [Angelopoulos, 2008]. The ESA measures 
plasma over the energy range of a few eV up to 30 keV for 
electrons and 25 keV for ions. The SST measures the distri- 
bution functions of superthermal particles within the energy 
range from 25 keV to 6 MeV for electrons and 900 keV for 
ions. Because the measured quasi-static electric field com- 
ponent along the spin axis has large uncertainty due to the 
short boom along the spin axis we use EB = to determine 
the electric field used to calculate the Ex B/B velocity. 
The coordinate systems used in this paper include geocentric 
solar magnetospheric (GSM) coordinates and field-aligned 
coordinates (FAC). The FAC is defined in the following way: 
The positive Z axis points in the direction of the magnetic 
field at the spacecraft's location. The positive X axis lies 
in the plane of the magnetic field line passing through the 
spacecraft's location, perpendicular to the Z axis, and points 
inwards (towards the inside of the field line). The positive Y 
axis completes the orthogonal right-handed system. Wave- 
forms are analyzed dynamically in time and frequency using 
Morlet wavelet transform [Torrence and Compo, 1998]. The 
wave vector is determined using Minimum Variance Analysis 
(MVA) [Khrabrov and Sonnerup, 1998] on bandpass filtered 
three magnetic field components. 



3. Observations 

Figure 1 shows a 7-min interval of the field and plasma 
observations made by probe E of the THEMIS mission on 
August 27, 2009. The boundary normal direction (deter- 
mined from MVA on the ambient magnetic field) was [0.99, 
0.015, -0.12] in GSM coordinates and almost identical to the 
GSM-X direction, consistent with the spacecraft being near 
the sub-solar point (indicated by the position parameters at 
the bottom of Figure 1). The spacecraft travels from the 
outer magnetosphere (SP) through the magnetopause (MP, 
indicated by two light green shaded bands) into the magne- 
tosheath (SH). The purple shaded band shows an ~12-sec 
interval of magnetic burst data capture. 

The magnetopause crossing can be seen in the change in 
Bz from positive to negative in Panel A of Figure 1. The 
differential energy flux of electrons in Panel K shows that in 



the magnetosphere, where Bz was positive, high-energy elec- 
trons were encountered; while in the magnetosheath, where 
Bz was negative, lower-energy electrons were measured. The 
spacecraft passed from the lower plasma density mangeto- 
sphere to the higher density (factor of 100) magnetosheath 
via a region of mixed magnetosheath/magnetospheric plas- 
mas comprising the low latitude boundary layer as shown 
in Panel H. This observation is evidence for the transport 
of solar wind plasma across the magnetopause. The pres- 
ence of accelerated plasma flow is seen through the mag- 
netopause current sheet as shown in Panels D and E. The 
spacecraft crossed the magnetopause south of the separator, 
as suggested by the negative GSM-Z component of ion flow 
velocity and the result of Walen test [Sonnerup et al., 1981]. 
Based on the prediction of the magnetopause reconnection 
model being a rotational discontinuity [Paschmann et al., 
1979], the Walen test [Sonnerup et al, 1981] states that 
the observed flow velocity change between a point in the 
magnetopause and a reference point in the adjacent mag- 
netosheath equals the predicted modified Alfven velocity 
change. The angle deviations between the observations and 
the prediction are almost 180° for this event. This anti- 
parallel relation indicates that the spacecraft crossed south 
of the separator [Sonnerup et al., 1981]. The encounter with 
the magnetopause current sheet is associated with fast ion 
jetting consistent with the Walen relation and fast electron 
flows, indicating that reconnection is occurring. Magnetic 
reconnection is generally considered to be the primary mech- 
anism through which transport of plasma and energy across 
the magnetopause occurs. 

3.1. Identification of EDR 

Enhanced wave activity can be seen associated with the 
magnetopause crossing from Panels A (ambient magnetic 
field), B (burst magnetic field) and C (electric fluctuations). 
We note that the electric fluctuations maximize during the 
magnetic burst interval. During this interval, electron flow 
speed (Panel E), anisotropy, agyrotropy and Mach number 
(Panel J) also maximize. These enhanced amplitudes are 
coincident with abrupt increases in electron parallel tem- 
perature T e || shown in Panel I suggesting strong wave heat- 
ing. As will be discussed in more detail in next section, the 
observed intense waves may provide the observed electron 
heating. All these features, along with the fact that the elec- 
tron perpendicular flow velocity is not consistent with the 
E x B/B 2 velocity during the magnetic burst interval (Panel 
G), provide evidence for the detection of an EDR. Another 
feature is the density depletion shown in Panel H (purple 
shaded band) which might indicate the center of the EDR. 
These signatures are consistent with the simulation and ob- 
servations of Scudder et al. [2012]. Panels L and M show the 
electron pitch angle spectra for lower engery electrons mea- 
sured by ESA and higher energy electrons measured by SST, 
respectively. The magnetic field (Panel A) and electric field 
(Panel C) fluctuations enhance in the magnetopause bound- 
ary layer and in the magnetosheath with field-aligned and 
counter-streaming lower energy electrons as shown in Panel 
L. However, the higher energy electron pitch angle (Panel 
M) enhances around 90° in the boundary layer and in the 
magnetosheath, especially during the purple shaded mag- 
netic burst interval. Distinct from the electron distributions 
in the boundary layer and in the magnetosheath, the higher 
energy electron pitch angle enhances at 0° and 180° and the 
lower energy electrons are more isotropic near the current 
sheet center around 15:36:00 UT and 15:37:00 UT. 
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3.2. Observations of waves 

Figure 2 shows an example of the identified WH waves at 
the time indicated by a black vertical line in Panel B of Fig- 
ure 1. It can be seen from Panels A and B of Figure 2 that 
the waves have frequencies from 0.1 to 0.6 / ce (electron cy- 
clotron frequency) with amplitudes up to 3 nT (peak-peak). 
As can be seen from Panel C, the wave Poynting flux is 
mostly anti-parallel to the ambient magnetic field. The wave 
vector (k = [-0. 23, -0.09, 0.97] in FAC), determined from 
MVA on bandpass filtered wave magnetic field, is nearly 
along the background magnetic field. The wave propaga- 
tion angle with respect to the ambient magnetic field 8 kB 
is determined as ~ 166° since the Poynting flux is mostly 
anti-parallel to the magnetic field. An expanded view of the 
WH waves can be seen in Panels D and E which respectively 
show the filtered (200-2000 Hz) waveforms of the burst mag- 
netic and electric field data over the time interval indicated 
by the purple bar in Panel A. Panel F shows that the WH 
waves are circularly right-handed polarized with respect to 
the ambient magnetic field as expected. The electron dis- 
tribution functions observed at times close to and/or con- 
currently with the WH waves are shown in Panels G and 
H. The lower energy electrons (~100 eV) shown in Panel G 
have anisotropic distributions with T e ±/T e « < 1. However, 
the energetic electrons (~30 keV) shown in Panel H have 
anisotropic distributions with a larger population moving 
in perpendicular direction. Broad-banded emissions with 
strong electric (~10 mV/m) and magnetic (~40 nT) field 
fluctuations below ion cyclotron frequency (~1 Hz) and elec- 
tric (~30 mV/m) and magnetic (~20 nT) field fluctuations 
below LH frequency (~30 Hz) are also detected during this 
magnetopause crossing (not shown). These intense wave 
emissions may provide the observed electron heating associ- 
ated with the magnetopause crossing. 

Figure 3 shows examples of electrostatic waves at a time 
preceding the waves in Figure 2 by 0.05 seconds. WH waves 
at a frequency of ~ 0.3/ ce are observed in the magnetic fluc- 
tuations in Panel A from 15:35:33.650 UT to 15:35:33.700 
UT. This time interval overlaps with that of the ESWs, sug- 
gesting a possible coupling of WH waves and ESWs. ESWs 
(up to 30 mV/m) indicated by the magenta arrows in Panel 
B have a broad spectrum which extends from 200 Hz to 
3000 Hz shown in Panel C. The high-frequency electrostatic 
waves (up to 35 mV/m) labeled by light blue arrows in 
Panel B have power that peaks at f ce , which can be seen 
from both the wavelet power spectrum in Panel C and the 
Fourier power in Panel G. During this time interval, there is 
no power in the wave magnetic field at f ce (not shown). This 
wave mode is linearly polarized, as shown in the hodograms 
in Panels D, E, and F in FAC with an interval indicated 
by the light blue arrows below Panel C. Occasionally these 
high-frequency emissions are seen associated with clear har- 
monics, possibly suggesting electrostatic electron cyclotron 
(EEC) waves. 



4. Discussion and Conclusions 

We have presented the first observations of intense waves 
in the EDR in a sub-solar magnetopause reconnection re- 
gion. The identification of the EDR in this event is based on 
the occurrence of signatures of strong electron heating, large 
electron thermal anisotropy, agyrotropy and Mach number 
and electron velocity not consistent with the E x B velocity, 
consistent with the simulation and observations of Scudder 
et al. [2012]. The lower energy electrons (~100 eV) with 



anisotropic distributions of T e ±/T e « < 1 within the EDR 
may have been heated by the observed waves with frequency 
below the LH frequency, consistent with the suggestion that 
LH waves lead to electron heating in the parallel direction 
[Cairns and McMillan, 2005]. 

We identified intense WH waves inside the EDR. This is 
inconsistent with reported simulation results that indicated 
WH waves are only driven downstream of an EDR [Fuji- 
moto and Sydora, 2008]. The WH waves seen by THEMIS 
in the EDR propagate almost anti-parallel to the ambient 
magnetic field and the Poynting flux indicates that the WH 
waves propagate away from the center of the 'X-line' along 
magnetic field lines. The observed electron temperature 
anisotropy of T e ±/T e u > 1 for energies above 20keV may 
be the source of free energy for the generation of the WH 
waves. The energetic electron anisotropy may be produced 
by adiabatic heating in the perpendicular direction as the 
locally intensified magnetic field can accelerate electrons in 
the perpendicular direction [Fujimoto and Sydora, 2008]. 
On the field lines directly connected to the EDR, the en- 
ergetic electron anisotropy may also be due to high energy 
field-aligned electrons (accelerated by the reconnection pro- 
cess) being lost to the magnetosheath [Stenberg et al., 2005]. 
WH waves can scatter the electrons in pitch-angle distri- 
bution and relax the temperature anisotropy. Studies of 
large amplitude whistlers in the inner magnetosphere have 
provided evidence for rapid scattering and/or energization 
[Cattell et al., 2008]. WH waves may play a significant role 
in the microphysics of reconnection through the enabling of 
a current sheet instability, the decoupling of electrons, the 
acceleration and heating of particles, and the transport of 
energy away from the reconnection region. 

A possible coupling of EEC waves and ESWs with WH 
waves is often seen during magnetopause reconnection. The 
growth of the electrostatic waves may reduce the electron 
temperature anisotropy and reduce the growth rate of WH 
waves. The physics of wave coupling process is important to 
understand the effect of wave-wave interactions on the re- 
connection process and will be investigated in a future study. 
This study provides further evidence that the plasma waves 
can play a significant role in the microphysics of magnetic 
reconnection at the Earth's magnetopause. 
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Figure 1. A reconnection event at the sub-solar mag- 
netopause observed by THEMIS-E on August 27, 2009. 
(A): 4 samples/s magnetic field data in GSM. (B): per- 
pendicular X component of the burst magnetic field at 
8192 samples/s in FAC. (C): perpendicular Y component 
of the electric field at 128 samples/s in FAC. (D) and (E): 
ion and electron bulk flow velocity in GSM, respectively. 
(F) and (G): comparisons of GSM-Z component of the 
E x B /B 2 velocity with the GSM-Z component of ion (F) 
and electron (G) perpendicular flow velocity with respect 
to the ambient magnetic field, respectively. (H): ion and 
electron densities. (I): electron temperatures. (J): elec- 
tron agyrotropy and Mach number (scale to the left) and 
temperature anisotropy (scale to the right). (K): differ- 
ential energy flux for electrons measured by ESA. (L) 
and (M): electron pitch angle spectra for lower engery 
electrons measured by ESA and higher energy electrons 
measured by SST, respectively. 
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Figure 2. Example of WH waves within the EDR. (A), 
(B) and (C): a 0.53-sec interval in FAC of perpendicu- 
lar X component of the burst magnetic field, associated 
Wavelet power spectrum and whistler Poynting flux, re- 
spectively. (D) and (E): expanded views of the filtered 
whistler waveforms of the perpendicular X component of 
the burst magnetic field and the perpendicular Y compo- 
nent of the electric field over the time interval indicated 
by the purple bar in (A). (F): hodogram of the filtered 
burst magnetic field waveforms in FAC for the interval 
indicated by the green bar in (D). The black star and 
dot mark the beginning and ending of the wave field, re- 
spectively. (G) and (H): distribution functions of lower 
energy electrons (up to 20 keV) measured by ESA and 
higher energy electrons (20-700 keV) measured by SST 
observed at times close to the WH waves, respectively. 
The horizontal axis is parallel to the ambient magnetic 
field and the bulk velocity defines the plane. 
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Figure 3. Example of electrostatic waves within the 
EDR. (A), (B) and (C): a 0.1-sec interval of the per- 
pendicular X component of the burst magnetic field, the 
parallel Z component of the electric field waveform cap- 
ture and the associated Wavelet power spectrum of the 
parallel electric field in FAC, respectively. (D), (E) and 
(F): hodograms of the electric field waveforms with an in- 
terval indicated by the light blue arrows below (C). The 
definitions of the black star and dot on hodograms are the 
same as those in Figure 2. (G): Fourier wave power vs. 
frequency with the same time interval as the hodograms. 



